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I.  Introduction 

The  motion  of  a  satellite  through  the  terrestrial  ionosphere  creates  dis¬ 
turbances  whose  spatial  extend  and  temporal  characteristics  are  not  yet  quanti¬ 
tatively  known.  Neither  are  the  processes  (e.g.  particle  acceleration,  insta¬ 
bilities,  wave  excitation)  acting  in  the  wake  zone  behind  the  satellite  under¬ 
stood  in  terms  of  varying  plasma  flow  and  body  parameters. 

While  there  is  an  awareness  that  such  disturbances  may  have  an  impact  on 
the  reliability  of  low-energy  particle  measurements,  we  find  that  most  published 
data  (from  in-situ  observations)  relating  to  such  disturbances  are  by-products 
of  geophysical  measurements  and  that  no  satellite  mission  had  such  a  study  as  a 
major  experimental  objective.  Hence,  the  available  information  is  fragmentary 
And  does  not  always  lend  itself  to  physically  meaningful  parametric  investigation 
and/or  theory-experiment  comparisons. 

Theoretical  investigations  of  the  effects  which  occur  in  the  near  and  far 
vicinity  of  spacecraft  moving  in  planetary  ionospheres  require  the  self-consistent 
solution  of  the  Boltzmann-Poisson  or  the  Vlasov-Poisson  (for  a  collisionless 
plasma)  equations  under  boundary  conditions  specified  by  both  spacecraft  and 
plasma  properties.  This  is  not  an  easy  problem.  Therefore,  in  order  to  make 
the  mathematical  treatment  more  tractable,  simplifying  assumptions  are  used, 
whose  validity  and  physical  significance  in  determining  the  overall  flow  field 
around  the  spacecraft  are  not  yet  clear.  A  discussion  of  various  aspects  of  the 
problem  including  those  specific  cases  which  lend  themselves  to  semi-analytical 
treatment  is  given  in  Al’pert,  1976.  A  discussion  of  theoretical  work  dealing 
with  various  numerical  methods  and  approaches  is  given  by  Parker  [1976,  1977). 
Briefly,  among  the  deficiencies  of  available  theoretical  models  is  the  use  of 
linearization  procedures  which  require  the  plasma  potential  of  the  disturbed 
regions  around  the  spacecraft  to  be  small  and  solutions  are  sought  at  large 
distances  from  the  spacecraft  (where  the  geometry  and  size  of  the  body  and 


^  V*—« 


sometimes  its  potential  are  less  significant).  In  other  models  the  positive 
ions  are  treated  as  neutral  particles  (i.e.,  the  "neutral  approximation")  thus, 
ignoring  the  effect  of  electric  fields  on  the  ion  trajectories.  Another  major 
deficiency  of  present  day  theories  is  that  the  non-stat ionary  case  has  not  been 
realistically  addressed.  For  discussions  on  various  aspects  of  the  problem, 
see  Al'pert  [1976];  Liu  [1975a];  Gurevich  and  Pitaevskii  11975] ;  Gurevich  and 
Dimant  [1975];  Gurevich  et  al.  [1973,  1968,  1966].  In  other  words,  despite  the 
fact  that  wave  excitation  and  plasma  instability  in  the  vicinity  of  a  moving 
body  seem  to  exist,  [Samir  and  Willmore,  1965;  Samir  and  Wrenn,  1969;  Gurevich 
et  al.,  1973;  Liu,  1975b;  Al'pert,  1976],  no  detailed  studies  on  this  basic 
aspect  of  the  problem  have  been  performed  .  The  latter  is  surprising  in  view  of 
the  potential  application  of  this  area  of  research  to  Astrophysics. 

Experimental  tests  of  assumptions  used  in  various  wake  models  are  practically 
complicated  since  most  of  the  available  in-situ  measurements  were  made  by  flush- 
mounted  probes,  hence  the  measurements  are  confined  to  the  very  near  vicinity  of 
the  bodies.  The  results  presented  in  this  paper  are  subject  to  the  same  limi¬ 
tation.  Since  the  self-consistent  numerical  solution  to  the  problem  of  body- 
plasma  interaction  (for  realistic  situations  in  space)  is  difficult,  it  seems 
that  parametric  and  theory-experiment  investigation,  limited  and  imperfect  as 
they  may  be,  can  contribute  significantly  to  a  better  physical  understanding  of 
the  interactions.  The  present  investigation  is  aimed  at  presenting  and  dis¬ 
cussing  some  experimental  results  which  are  useful  to  some  scientific  aspects 
of  satellite-ionosphere  interactions,  as  well  as  to  technological/application 
aspects.  In  particular  to  the  planning  of  future  experiments  of  body-plasma 
electrodynamic  interactions  in  a  supersonic  and  subAlfvenic  flow  regime  to  be 
conducted  on  board  the  Shuttle/Spacelab  and  via  the  utilization  of  subsatellites 
and  tethered  bodies.  More  specifically,  it  is  expected  that  investigations  of 


the  kind  presented  here  be  useful:  (1)  In  the  planning  of  instrument  location 
on  ejectable  ensembles  of  probes  and  on  the  orbiter  itself  in  future  shuttle 
missions,  (2)  in  predicting  the  disturbances  expected  to  be  created  by  deployed 
tethered-balloons  and  other  boom-mounted  bodies  in  future  shuttle  missions  and 
(3)  in  testing  theoretical  wake  models.. 


One  of  the  quantities  that  can  serve  as  a  measure  of  the  degree  of  distur- 

I+(wake) 

batnce  produced  by  the  spacecraft  motion  is  the  ratio  a  =ly - ].  Unfor- 

+(ambient) 

tunately,  there  is  little  systematic  experimental  information  regarding  the 

Ro  e4>s 

dependence  of  a  on  basic  plasma  parameters  such  as:  R^  =  Cj^)  ;  <|>^  =  (j^— ] )  ; 

S  =  s  .)  and  (T  /T.)  where:  R  =  satellite  radius,  =  Debye  length,  <t>  = 

12  KT  _  e  +  O  D  1  6  s 


satellite  potential,  =  electron  temperature,  =  ion  temperature,  V  =  satellite 
velocity,  M+  =  ionic  mass,  K  =  Boltzmann's  constant,  e  =  electronic  charge. 

It  is  the  major  objective  of  the  present  investigation  to  contribute  to  our 


knowledge  and  understanding  of  the  variation  of  a  with  some  plasma  and  body 


parameters. 


II.  Experimental  -  The  Instruments  Used 

In  order  to  assess  the  relative  Influence  of  the  parameters  S,  and  $ 
on  the  amount  of  ion  depletion  in  the  wake  of  the  S3-2  satellite,  measurements 
of  ion  current  (I+)  electron  temperature  (T£)  ,  electron  density  (N^)  and  values 
of  plasma  potential  ($g)  from  the  planar  multigrid  probes  and  from  a  spherical 
gridded  probe  were  used.  The  planar  ion  probes  (or  sensors)  were  flush-mounted 
on  the  surface  of  the  satellite  and  the  spherical  electron  probe  was  mounted  at 
the  end  of  a  boom  1.3  m  long  and  parallel  with  the  satellite  spin  axis.  The 
ion  probes  were  actually  contained  in  two  packages  placed  180°  apart  and  with 
centerlines  in  the  spin  plane  of  the  satellite.  Each  package  contains  four 
independent  planar  ion  sensors.  Two  of  these  sensors  view  in  the  spin  plane  of 
the  satellite,  and  two  look  at  40°  on  either  side  of  the  spin  plane.  The 
principal  objective  of  the  sensor  geometry  was  to  determine  direction  of  bulk 
ion  flow  relative  to  the  spacecraft  when  the  sensor  array  is  close  to  the  'ram' 
direction.  However  the  'ram'  current  flowing  to  the  sensors,  together  with 
current  measured  at  other  points  in  the  satellite  rotation  may  be  used  to  help 
build  a  picture  of  the  positive  ion  distribution  around  the  satellite.  The 
location  of  the  planar  ion-probes  on  the  satellite  as  well  as  the  location  of 
the  electron  probe  is  shown  schematically  in  Figure  1  (a  and  b) .  Details  of  the 
ion-probe  internal  configuration,  electronics,  probe  materials  and  a  discussion 
on  probe  performance  is  given  in  Wildman  [1977].  A  detailed  discussion  on  the 
instrument  package  on  the  S3-2  satellite  is  given  in  Wildman  [1976]  and  Lai 
et  al.  [1976],  and  will  not  be  further  discussed  here. 
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If  a  *  f(R^)  at  R^  *  40  (from  Figure  4)  depicts  the  behavior  of  a  due 
t0  an<^  not  ^ue  to  an  1  to  (H  ]  transition  (which  happened  to  occur 
at  the  same  altitude)  than  it  is  possible  that  the  behavior  of  a  ■  f(R^)  for 
Rp  >_  40  signifies  the  behavior  due  to  the  ’large-body'  case  [e.g.  Al'pert,  1965; 
Call,  1969;  Liu,  1969;  Parker,  1976;  Al'pert,  1976]  where  a  significant  drop 
in  a  is  expected  to  occur  very  close  to  the  satellite's  surface.  At  this  point 
it  is  not  unreasonable  to  interpret  the  variation  in  a  including  the  abnpt  drop 
in  terms  of  o  =  f(R^)  mainly.  However  the  possibility  that  this  interpretation 
is  not  unique  exists,  e.g.  there  could  have  been  an  [0+]  to  [H+]  transition  at 
H  =  500  km.  The  measurements  from  the  S3-2  satellite  could  not  have  revealed 
such  a  tra  Ltion,  since  no  ionic  composition  measurements  were  available  to  us. 
Figure  5  shows  the  variation  of  a  =  f(Rp)  for  25  <_  R^  £  56.  As  could  have  been 
expected  [e.g.  Samir  et  al.,  1976]  a  attains  an  exponentially  decreasing  value 
for  increasing  values  of  R^. 

Using  S3-2  data  we  have  examined  a  “  f(<)>N)  for  |<}>N[  in  the  range  10  <  |  [  < 

18.  The  results  are  shown  in  Figure  6.  As  seen,  no  obvious  correlation  between 
the  abrupt  change  in  a  =  f(H)  can  be  attached  to  the  a  =  f(<|>N)  variation.  In 
fact,  as  seen  from  this  figure,  significantly  different  values  of  a  are  obtained 
for  the  same  value  of  which  may  imply  a  behavior  due  to  the  influence  of 
another  parameter.  Examination  of  the  results  shown  in  Figure  5  show  that  lower 
values  of  a  correspond  to  higher  values  of  and  higher  values  of  a  correspond 
to  lower  values  of  R^  for  similar  values  of  <t^.  This  may  suggest  that  very  close 
to  the  satellite  surface,  R^  dominates  a  more  than  does 

As  mentioned  earlier  we  have  attempted  to  assess  the  relative  importance  of 
SAV  in  determining  the  ion  distribution  in  the  very  near  wake  and  examine  the 
importance  of  SAV  vis-a-vis  that  of  R^.  In  carrying  out  a  traditional  analysis 
and  using  data  for  [Mf (average) ]  from  the  literature  [e.g.  Banks  and  Kockarts, 
1973;  Al'pert,  1973]  we  could  find  no  obvious  abrupt  change  in  S  ^  that  could  be 


correlated  with  a  =  f(H)  at  about  500  km.  The  change  in  is  being  rather 
gradual. 

Since  no  local  ion  composition  measurements  from  the  S3-2  satellite  were 
available  we  have  attempted  to  use  ion  measurements  from  another  satellite  with 
orbits  close  enough  to  the  ones  we  used  from  the  S3-2  satellite.  It  appeared 
that  some  orbits  of  the  Atmosphere  Explorer  E  (AE-E)  satellite  could  be  con¬ 
sidered  for  this  purpose.  Therefore,  we  have  examined  a  fairly  large  sample  of 
ion  composition  measurements  from  the  Bennett  ion  mass  spectrometer  (BIMS) 
mounted  on  board  the  AE-E  satellite. 

The  mode  of  operation,  technique  and  method  used  in  the  BIMS  experiment  are 
given  in  detail  elsewhere  [Brinton  et  al . ,  1973]  and  will  not  be  discussed  here. 
Figure  7  shows  an  example  of  the  final  format  of  the  measurements  we  obtained  from 
the  BIMS  experiment.  In  this  figure  the  detailed  variation  of  [0+] ,  [H+] ,  [N+] , 

[He+]  densities  with  altitude  are  shown.  From  such  information  the  transition 

4-4*  4- 

altitudes  of  [0  ]  to  [H  ]  is  readily  obtained.  In  fact  the  results  of  [N(0  ), 

N(H+),  N(N+),  N(H  +) ]  =  f(H)  shown  in  Figure  7  are  from  the  closest  orbits  we 
e 

could  match  to  the  S3-2  orbit  whose  a  =  f(H)  results  are  given  in  Figure  3.  It 
should  be  realized  that  matching  'similar*  orbits  for  the  S3-2  and  the  AE-E 
satellite  is  subject  to  limitations.  This  is  due  to  the  fact  that  in  almost  no 
case  could  we  reach  a  situation  where  good  coverage  of  'overlapping'  orbits  was 
possible.  For  example,  latitude  could  be  matched  to  within  (20-25°),  longitudes 
had  a  wide  diversity  and  local  time  could  be  matched  to  within  (3-5)  hours 
difference  or  more.  Figure  7  depicts  as  reasonable  a  match  as  possible  between 
AE-E  ion  composition  results  and  those  shown  in  Figure  3  for  S3-2  results.  As 
seen  from  Figure  7  a  difference  of  H  =  150  km  exist  between  the  [0+]  to  [H+] 
transition  altitudes  for  the  two  orbits  shown. 
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It  was  interesting  to  determine  the  ratio  of  normalized  densities  [-  — ) 

o 

where:  *=  is  the  ambient  density  and  90°  6  <_  180°  at  the  [0+]  to  [H+] 

transition  altitude  where  N(0+)  =  N(H+)  and  then  see  where  does  the  value  of  these 
ratios  fall  on  the  a  =  f(H)  plots  e.g.  Figure  3.  For  [N(0)/Nq)]  we  have  used 
expression  (1)  given  below. 

N(6)  _  n(Q+  ^1+erf [S(0+)  -K-cosO]  ^  +  ^1+erf  [S(H~l')K*cos8]  ^  ^ 

No  1+erf [S(0+)  K]  l+erf[S(H+)  K] 


In  this  expression  [after:  Gurevich  et  al.,  1970;  Al'pert,  1976]  K  is  a  constant 

(:  0.6),  6  =  angle  of  attack,  N(0)  =  density  at  the  angle  0 ,  =  ambient  density, 

n(0+)  and  n(H+)  =  relative  concentration  of  the  [0+]  and  [H+]  ions  in  the  plasma 
+  + 

and  S(0  )  and  S(H  )  =  corresponding  ionic  Mach  numbers  considered  separately. 

In  fact,  in  obtaining  expression  (1)  the  two  ion  species  are  considered  non¬ 
interacting  and  the  influence  of  the  electric  field  upon  ion  trajectories  small. 
Somewhat  more  detailed  discussions  on  the  latter  are  given  in  the  Appendix. 

In  computing  ~ ~]  we  have  used  the  ion  measurements  from  the  BIMS  experi- 

o 

ment  [Brinton  et  al.,  1973]  on  the  AE-E  satellite  with  electron  temperature 

measurements  from  the  spherical  probe  on  the  S3-2  satellite.  We  found  that  the 

values  of  having  ratios  similar  to  a  ‘occur  sometimes  close  to  the  altitude 

o 

where  the  abrupt  drop  in  a  occurs  and  sometimes  at  altitudes  faraway  (by  hundreds) 

n(g)  i(g) 

of  kms)  from  this  altitude.  Indeed  by  comparing  [  —  - — ]  values  with  [— -  --]  values 

o  o 

we  impose  a  limitation  on  the  study,  and  this  subject  was  discussed  in  detail  in 
Samir  et  al.  [1979a]  and  Samir  et  al.  [1980].  At  this  point  and  despite  our 


extensive  effort  we  cannot  claim  that  we  have  separated  between  effects  due  to 
Rp  and  those  possibly  due  to  the  transition  of  [0+]  to  [H+] .  On  the  other  hand 
we  have  observed  that  rather  than  using  the  parameter  S ^  it  is  more  useful  and 
realistic  to  use  actual  [M+(specif ic) ]  values  together  with  specific  values. 
This  is  in  accord  with  Samir  et  al.  [1979b]  and  with  the  discussion  given  in 


Al'pert  [1976]  and  Gurevich  and  Pitaevskii  [1975]  regarding  the  impoitance  of 

percentage  [H+]  in  determining  the  ion  and  electron  distribution  in  the  wake. 

From  other  studies  [e.g.  Samir  et  al . ,  1979a;  Samir  et  al. ,  1980]  we  know  that 

u.  =  f(R_)  and  that  a  =  f  (^  ^ ,  T  )  [e.g.  Samir  et  al .  ,  1979b]  hence  a  dependence 
o  N (H+)  e 

on  both  R^  and  S(0  )  and  S(H+)  does  exist.  On  the  other  hand  the  use  of  'average' 
parameters  (e.g.  S^)  should  be  exercised  with  care  and  reservation  depending  on 
the  main  objective  of  the  study.  For  example  the  S  ^  is  more  influenced  by  the 
[0+]  ions  than  by  [H+]  ions  having  the  same  density.  The  latter  will  have  a 
.  imilar  influence  on  S.y  as  will  [N+]  ions  having  a  density  which  is  one  order 
of  magnitude  less  than  the  [H+]  density. 

It  is  possible  that  in  the  very  near  wake  the  parameter  =  f(T  ,  N^) 
competes  with  the  specific  ionic  Mach  number  (S)  rather  than  with  [M+(average) ] 
or  sAy  That  a  is  indeed  a  function  of  has  been  shown  in  Samir  et  al .  [1979a], 
Samir  [1980]  and  the  relationship  a  =  f(Rp)  was  shown  to  be  exponential  for  R^  >_ 
(40-50).  Most  of  the  results  dealt  with  in  the  above  papers  were  confined  to 
altitudes  which  are  below  the  [0+]  to  [H+]  transition  with  the  possible  exception 
of  case  if 4  [  Table  2  in  Samir  et  al.,  1979a].  And  indeed  in  analyzing  the  above 
case  the  possible  contribution  of  [N(H+) /N(0+) ]  was  mentioned. 

From  the  present  study  it  is  obvious  that  more  parametric  investigations 
are  needed  in  order  to  clarify  in  a  more  quantitative  manner  the  parametric 
interplay  which  determines  the  ion  distribution  around  the  satellite.  Such  studies 
are  also  directly  relevant  to  the  determination  of  the  surface  of  the  potential 
on  bodies  in  space  [e.g.  Whipple,  1977]. 
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IV.  Suranary 

The  variation  of  Ion  current  depletion  in  the  wake  of  the  S3-2  satellite 
with  altitude  in  the  range  300  km  to  1100  km  was  quantitatively  determined.  An 
abrupt  drop  in  ion  distribution  at  about  500  km  can  be  correlated  with  either 
a  significant  change  in  and/or  with  a  transition  between  [0+]  to  [H+]  or  with 
a  combination  of  both.  The  relative  importance  of  the  parameter  versus 
in  determining  the  ion  distribution  in  the  wake  was  examined  quantitatively  in 
the  range:  10  <  I  <{•  [  <18. 

Parametric  studies  such  as  attempted  here  should  be  continued  on  board  the 
Spacelab/Orbiter  if  we  are  to  uniquely  resolve  the  intricacies  of  the  parametric 
interplay  in  body-plasma  interactions,  in  supersonic  and  subAlvenic  flow  regimes. 


Acknowledgement 

This  work  was  supported  by  the  U.S.  Air  Force  contract  F19628-79-C-0060 . 
We  thank  Dr.  Irving  Michael  from  the  Air  Force  Geophysics  Laboratory  and 
Mr.  John  Metzger  from  the  Space  Physics  Research  Laboratory,  University  of 
Michigan  for  their  help  in  handling  S3-2  and  AE-E  data. 


i- 


References 


Al'pert,  Y.  L.,  A.  V.  Gurevich,  and  L.  P.  Pitaevskii,  Space  physics  with 
artificial  satellites,  Monography,  Consultant  Bureau,  1965. 

Al'pert,  Y.  L.,  Radio  wave  propagation  and  the  ionosphere.  The  Ionosphere, 

Vol.  1,  Consultant  Bureau,  1973. 

Al'pert,  Y.  L.,  Wave-like  phenomena  in  the  near  earth  plasma  and  interactions 
with  man-made  bodies,  Handbuch  der  Physik  (ed.  S.  Flugge)  Geophysics  III, 
Part  V,  217,  1976. 

Banks,  P.  M. ,  and  G.  Kockarts,  Aeronomy ,  Academic  Press,  1973. 

Banks,  P.  M. ,  R.  W.  Schunk,  and  W.  J.  Raitt,  The  topside  ionosphere:  a  region 
of  dynamic  transition.  Annual  review  of  earth  and  planetary  sciences,  _4, 
381,  1976. 

Brinton,  H.  E. ,  H.  G.  Mayr,  R.  A.  Pickett,  and  H.  A.  Taylor,  Jr.,  The  effect 

of  atmospheric  wind  on  the  0+-H+  transition  level.  Space  Res.  X,  652,  1970. 

Brinton,  H.  C. ,  L.  R.  Scott,  M.  W.  Pharo,  and  J.  T.  C.  Coulson,  The  Bennett  ion- 
mass  spectrometer  on  Atmosphere  Explorer  C  and  E,  Radio  Sci. ,  8,  323,  1973. 

Call,  S.  M. ,  The  interaction  of  a  satellite  with  the  ionosphere.  Report  No.  46, 
Plasma  Laboratory,  Columbia  Univ. ,  1969. 

Gurevich,  A.  V.  ,  L.  V.  Parilskava,  and  L.  P.  Pitaevskii,  Self-similar  motion 
of  rarefied  plasma  I,  Soviet  Physics/JETP,  22(2) ,  449,  1966. 

Gurevich,  A.  V.,  L.  V.  Pariiskaya,  and  L.  P.  Pitaevskii,  Self -similar  motion 
of  a  low  density  plasma  II,  Soviet  Physics/JETP,  27(3) ,  476,  1968. 

Gurevich,  A.  V.,  L.  P.  Pitaevskii,  and  V.  V.  Smirnova,  Ionospheric  aerodynamics, 
Soviet  Physics  -  Uspekhi,  £9,  (1-2),  595,  1970. 

Gurevich,  A.  V.,  L.  V.  Pariiskaya,  and  L.  P.  Pitaevskii,  Ion  acceleration  upon 
expansion  of  a  rarefied  plasma,  Soviet  Physics/JETP,  36(2) ,  274,  1973. 

Gurevich,  A.  V.  ,  and  Ya.  S.  THmant,  Flow  of  a  rarefied  plasma  around  a  disc, 
Geomag.  and  Aeron.,  183,  1975. 


i  ■  J  ^  i  *  • 


Gurevich,  A.  V.,  and  L.  P.  Pitaevskii,  Non-linear  dynamics  of  a  rarefied  Ionized 


gas,  Prog.  Aerospace  Sci.,  16 ,  227,  1975. 

Lai,  S.  T.,  M.  Smiddy,  and  P.  J.  L.  Wildman,  Satellite  sensing  of  low  energy 

plasma  bulk  motion,  7th  Conference  on  Aerospace  and  Aeronautical  Meteorology 
and  Symposium  on  Remote  Sensing  from  Satellites,  November  16-19,  1976. 

Liu,  V.  C.,  Ionospheric  gas  dynamics  of  satellites  and  diagnostic  probes. 

Space  Sci.  Rev.,  _9,  423,  1969. 

Liu,  V.  C.,  On  ionospheric  aerodynamics,  Prog.  Aerospace  Sci.,  16,  273,  1975a. 

Liu,  V.  C.,  A  wave  model  of  near  wakes,  Geophys.  Res.  Lett.,  2,  485,  1975b. 

Parker,  L.  W.,  Computation  of  collisionless  steady-state  plasma  flow  past  a 
charged  disc,  Report,  NASA  CR-144159,  February  1976. 

Parker,  L.  W. ,  Calculation  of  sheath  and  wake  structure  about  a  pillbox-shaped 
spacecraft  in  a  flowing  plasma,  Proc.  Spacecraft  Charging  Technology 
Conference  (ed.  C.  P.  Pike  and  R.  R.  Lovell)  p.  331,  1977. 

Samir,  U.,  and  A.  P.  Willmore,  The  distribution  of  charged  particles  near  a 
moving  spacecraft,  Planet.  Space  Sci.,  13 ,  285,  1965. 

Samir,  U.,  and  G.  L.  Wrenn,  The  dependence  of  charge  and  potential  distribution 
around  a  spacecraft  on  ionic  composition.  Planet.  Space  Sci.,  17 ,  693,  1969. 

Samir,  U.,  E.  J.  Maier,  and  B.  E.  Troy,  Jr.,  The  angular  distribution  of  ion 
flux  around  an  ionospheric  satellite,  J.  Atmos.  Terrest.  Phys. ,  35 ,  513, 
1973. 

Samir,  U. ,  R.  Gordon,  L.  H.  Brace,  and  R.  F.  Theis,  The  near  wake  structure  of 
the  Atmosphere  Explorer  C  (AE-C)  satellite  -  a  parametric  investigation, 

J.  Geophys.  Res.,  84 ,  513,  1979a. 

Samir,  U.,  L.  H.  Brace,  and  H.  C.  Brinton,  About  the  influence  of  electron 

temperature  and  relative  ionic  composition  on  ion  depletion  in  the  wake 
of  the  AE-C  satellite,  Geophys.  Res.  Lett.,  101,  1979b. 


* 


Samir,  U. ,  Y.  J.  Kaufman,  L.  H.  Brace,  and  H.  C.  Brinton,  The  dependence  of  ion 
density  in  the  wake  of  the  AE-C  satellite  on  the  ratio  of  body  size  to 
Debye  length  in  an  [0+]  -  dominated  plasma,  J.  Geophys.  Res..  85,  1769, 

1980. 

Titheridge,  J.  E.,  Ion  transition  heights  from  topside  electron  density  profiles, 
Planet.  Space  Sci.,  2A,  229,  1976. 

Whipple,  E.  C.  Jr.,  Modeling  of  spacecraft  charging,  Proc .  of  Spacecraft  Charging 
Technology  (ed.  C.  P.  Pike  and  R.  R.  Lovell),  225,  1977. 

Wildman,  P.  J.  L. ,  Studies  of  low  energy  plasma  motion.  Results  and  a  new 
technique,  AFGL-TR-76-0168  Report/Air  Force  Geophysics  Laboratory  - 
Space  Science  Division,  1976. 

Wildman,  P.  J.  L.,  A  low  energy  ion  sensor  for  space  measurements  with  reduced 
photo-sensitivity,  Space  Sci  ■  Inst ■  ,  3,  363,  1977. 


I 


] 


DIAGRAM  OF  SATELLITE  AND  PROBE  LOCATION 
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Figure  1: 


Schematic  drawing  showing  the  location  of  the  ion  and  electron 
probes  on  the  S3~2  satellite. 
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An  example  showing  the  general  behavior  of  the  ion  current  with 
time  for  the  altitude  ranges:  (a)  H  -  109^-1056  km,  (b)  H  -  5 1 0-477  km, 
(c)  338-317  km,  for  equatorial  and  middle  latitudes,  at  local  time 
22. 5h  ±  0.5h 


jt  rwm* 


**  •  f(M)  EOR  SEVERAL  IONIC  CONSTITUENTS 
<*E-£  SATELLHE/BiMS  EXPERIMENT) 


?*Vn  *.  . 

*>  i^J3.  • 

♦  0^°f  . 


......  .*> 


APPENDIX:  ABOUT  THE  EXPRESSION  USED  IN  THE  PAPER* 

0 


I .  About  Bibliography 

In  the  paper  (Gurevich  and  Pitaevskii,  1975)  a  detailed  discussion  of 
non-linear  dynamics  regarding  rarefied  ionized  gases  is  given  in  detail. 

In  Gurevich  et  al.  (1970)  and  Al’pert  (1976)  detailed  discussions  are  given 
regarding  the  various  aspects  of  the  interaction  between  a  rarefied  plasma 
and  a  supersonically  moving  body.  The  latter  two  references  extend  studies 
performed  by  Gurevich,  Pitaevskii,  Al'pert.and  colleagues  during  the  past  two 
decades  and  are  in  fact  a  continuation  to  Al'pert  et  al.  (1965),  which 
summarized  earlier  work  by  the  above  Russian  group. 

It  is  somewhat  surprising  that  while  an  extensive  theoretical  effort  took 
place  in  Russia  regarding  body-plasma  interaction  less  intensive  and  meaningful 
work  was  performed  in  the  U.S.A.  It  appears  justified  to  say  that  many 
theoretical  modeling  as  well  as  experimental  work  regarding  'spacecraft 
charging’  would  have  benefited  significantly  if  more  attention  would  have  been 
paid  to  the  Russian  theoretical  work. 

In  the  following  sections  we  deal  with  some  aspects  of  body-plasma  inter¬ 
actions  based  on  Russian  work. 

II .  Assumptions  and  Basic  Equations 

A.  We  first  deal  with  one-dimensional  quasi-neutral  flows  in  the  absence  of 
a  magnetic  field. 

The  basic  mathematical  approach  relies  on  the  analogy  to  hydrodynamics, 
in  particular  to  the  self-similar  solutions  which  do  not  contain  any  charac¬ 
teristic  dimensions  in  the  initial  and  final  conditions.  Namely,  the  time  t 
and  the  coordinate  x  appear  in  the  solution  of  such  problems  only  in  some 
combinations  of  (x/t) .  It  is  argued  (e.g.  Gurevich  and  Pitaevskii,  1975)  that 

The  contents  of  this  Appendix  is  part  of  a  review  paper  by  U.  Samir  now  in 
preparat ion. 
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self -similar  solutions  in  hydrodynamics  describe  a  large  class  of  physical 
problems  among  which  is  the  study  of . the  wake  behind  a  body.  Therefore  it 
should  be  of  interest  to  examine  the  situation  for  a  rarefied  plasma  described 
by  the  collisionless  kinetic  equation  with  a  self-consistent  field. 

In  this  appendix  we  summarize  some  of  the  findings  which  are  relevant  to 
body-plasma  interactions  in  planetary  ionospheres  and  magnetospheres,  in 
particular,  the  validity  of  the  approach.  This  will  be  done  via  the  comparison 
with  in  situ  experimental  results. 

Considering  the  one  dimensional  case  and  assuming  that  at  an  initial  time 
the  plasma  occupies  the  half -space  x  <  0  and  then  expands  into  the  other  half 
of  the  space  (x  >  0)  tgfgn  the  plasma  is  described  by  the  kinetic  equation  for 
the  distribution  function  for  ions: 


where  M«  mass  of  the  ions,  $  •  the  electric  potential.  ^  Note:  A  similar 
equation  can  be  written  for  the  electron  distribution  function^, 
and  by  the  Poisson  equation: 


hi  =  A/O.-fc)  -  J-f dv 


- 

where  N  =  ion  number  density,  Ng  *  electron  density.  We  also  assume  that: 

A4-  ^  H  ) 


(3) 


,'A 


i.e.,  the  equilibrium  distribution  of  ions  is  described  by  the  Boltzmann 
formula  (3).  Note :  There  are  other  distributions  for  electrons  (e.g. 
Gurevich  and  Pitaevskii  (1975),  p,  258.  We  now  adopt  the  Russian  notation 
whereby  is  given  in  eV,  hence 


(3) 


For  a  quasi-neutral  plasma,  and  the  argument  is  that  this  is  the  case  for  the 
greater  part  of  the  region  behind  a  body  (i.e.,  of  the  wake)  we  have: 

N  -  K/e.  w 

and  it  follows  that 

€  < f  =  r-  -L 

(5) 

which  then  yields  for  (1)  the  equation: 


Introducing  non-dimensional  and  self-similar  quantities: 


We  obtain  the  equation  for  the  ion-distribution  function: 


!  Dg_  t)  ^ 
£)ti  t)x 


with : 


(8) 


(9) 


For  the  case  <<  T^,  the  electrons  play  a  small  part.  Hence  in  this 
case,  the  last  term  in  equation  (1)  can  be  neglected  and  equation  (1)  becomes: 


(10) 


And  equation  (10)  is  the  equation  describing  a  free  molecular  flow  of  neutral 
gas. 

This  case  is  thus  called  the  "neutral  approximation".  The  general 
solution  of  (10)  is: 

--  Zr)  (ID 

where  f  (x,  v)  is  the  initial  distribution  function, 
o 

For  the  case  T&  >>  T^  the  part  played  by  the  electric_  fieW  becomes  very 
important .  To  first  approximation,  the  thermal  motion  of  the  ions  can  be 
neglected  (T  /7 ^  -*■  0)  ,  hence  the  distribution  function  for  ions  can  be  repre¬ 
sented  by  a  £-function  of  velocity,  i.e.,  has  the  form: 


(12) 
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Substitution  of  (12)  in  (1)  yield: 

&  4 ( 

The  eq.  (13)  coincide  with  the  hydrodynamic  equations  for  an  ideal  isothermal 
gas.  Note:  In  Gurevich  and  Pitaevskii  (1975),  Appendix  III,  p.  261  a  dis¬ 
cussion  on  the  inclusion  of  the  ion  thermal  motion  is  given. 

In  summary :  The  basic  idea  in  the  Russian  work  discussed  above  is  the 
possibility  of  using  a  self-similar  approach  i.e.,  seeking  solutions  which 
depend  on  (x/t)  in  a  way  which  is  analogous  to  (in  some  cases)  ideal  fluids 
in  hydrodynamics.  Therefore,  non-dimensional  variables  (e.g.  t,  g,  u)  are 
introduced  and  cases  of  interest  to  space  plasma  physics  are  treated  in  an 
analogous  manner. 

B .  Discussion  of  the  self -similar  method 
The  basic  equations  are: 

(U- r)£L-l-^L^=c>  ( 

^  or  A  ju  D-t  °  (  (8_9. 

'f  =  =  JL  M  \ 

le. 

v\  — •  f  0  v  0 


n* 


(8-9) 


In  mathematical  terms  it  is  possible  to  discuss  the  boundary  conditions  of 


M-*v£-I2\>L2 


[ JL  Y))  = 


Note:  Eq.  (14)  is  essentially  the  same  as  6)  in  the  preceding  section. 


in  terms  of  a  decay  of  an  arbitrary  initial  discontinuity.  Then  the  boundary 
conditions  for  (14)  for  an  arbitrary  initial  discontinuity  are: 


■PH., 


-f'-f  (vO  *  4  ^ 

x>£) 


where:  f  and  f  are  arbitrary  initial  distribution  functions  of  the  ions  on 
1  2 

the  left  hand  side  and  on  the  right  hand  side,  of  the  discontinuity.  In  terms 

x  M  ^ 

of  the  non-dimensional  variable  (T)  defined  as:  *0 m  ]  (defined  in 

t 

eq.  (7)),  we  obtain  for  the  distribution  functions  and  g^: 

And  for  and  chosen  as  Maxwellian  functions  in  the  form: 

f'  (-  jpj]  C 


Note:  The  ions  in  region  I  are  considered  here  to  be  at  rest, 
Hence  for  and  we  obtain: 


§,  =  ["  ^ 
n  -  I JtV H  - 


vnere:  u 


s  3  .  R  _JZL 

2tJ  y  r  7-  -r„ 


Equations  (8-9)  determines  g(t,  u)  via  the  method  of  characteristics 
i.e.,  considering  the  curves  on  which  g  has  a  constant  value  (see  Gurevich 


and  Pitaevskii,  1975).  The  equation  of  the  characteristics  is: 


du  .\  l  fCt)  \  ,  r-s 


(IS) 


where  F(t3  is  a  non-dimensional  force  acting  on  the  ions  which  is  proportional 
to  the  strength  of  the  electric  field  in  the  plasma. 

In  dealing  with  the  character  of  n(T)  for  the  "neutral  approx imat ion"  case 
Gurevich  and  Pitaevskii  (1975)  show  that  for  the  boundary  conditions  (15)  we 
get: 


(20) 


The  extremium  of  the  function  nCO  is  given  by  the  condition:  [dn/dc]  «  0, 


namely : 

-  fat*)  =o 


Now,  depending  on  the  relations  between  No,  ,  T^,  T vq  a  distinction  should 
be  made  for  the  following  three  cases: 


CO  ~n 

OO 

On)  "E  -rii~T, 

These  cases  are  being  discussed  and  the  behavior  of  the  distribution  function 
g  *  f(u)  is  given  for  each  of  the  above  cases  (see  Fig.  1,  p.  231  in  Gurevich 
and  Pitaevskii,  1975).  Also,  the  behavior  of  the  non-dimensional  force  F(t) 
is  given  for  each  of  the  above  cases  (see  Fig.  2,  p.  231  in  the  above  reference) 
In  Summary:  It  is  found  that  for  case  (i)  we  have: 

('  1  ~  ^  Ul)  (23) 

For  case  (ii)  we  have: 

/f  ~  **) 

Note :  The  above  analysis  is  for  the  ’’neutral  approximation". 

However,  with  the  exception  of  the  case  T^/T^  ■*  “  qualitative  conclusions 
are  valid  not  only  for  the  "neutral  approximation:  situation. 

For  the  asymptotic  behavior  of  the  function  (^0  for  X  +  +  ®,  particles 
with  large  velocities  (u)  are  most  important.  And  in  all  cases  except  for  the 
case  of  the  expansion  of  a  plasma  into  a  vacuum  (see  curve  2  in  Fig.  2c)  the 
jump  of  the  electric  potential  |  $  ^  ‘  +  +(r!  appears  to  be  finite  hence 

the  distribution  of  fast  particles  does  not  vary  with  the  field. 

Therefore,  the  asymptotic  expressions  (23)  and  (24)  remain  valid  when  the 
electric  field  is  taken  into  account.  Note:  This  is  a  very  significant 
conclusion  that  may  be  used  for  pratical  cases  and  simplify  the  calculations. 
Furthermore,  Gurevich  and  Pitaevskii  (1975)  claim,  that  the  qualitative  con- 


'.•*  JL  ■ 


siderations  concerning  the  shape  of  n(t)  remain  valid.  And  even  for  the  case 

where  the  electric  field  can  change  the  asymptotic  value  of  the  coefficient 
,  Te  ? 

of  Iexp(-  —  u^) ]  the  monotonic  behavior  of  [n (Zj  ]  does  not  change. 


Ill .  Discussion  of  Some  Relevant  Specific  Cases 

fy  'Bf.  The  free  expansion  of  a  plasma  into  a  vacuum 

Here  as  in  the  discussion  in  Section  II  the  problem  is  dealt  with  in  a 

self -similar  manner  (for  a  collisionless  plasma) .  The  density  and  the  velocity 

distribution  of  the  ions  are  obtained  (Gurevich  et  al.,  1968).  The  main  finding 

is  that  in  the  course  of  filling  the  rarefied  half-space  (e.g.  the  wake  region 

behind  a  body)  a  part  of  the  ions  is  accelerated  by  the  action  of  the  resulting 

electric  field  up  to  velocities  of  the  order  of  the  thermal  velocity  of  the 

electrons.  At  the  same  time  the  effective  temperature  of  the  ions  drops  sharply. 

In  fact  the  drop  Is  such  that  the  effective  ion  temperature  is  many  times  smaller 

than  the  electron  temperature  (for  the  initial  situation  of  T  *  T  ). 

6  * 

Here  we  will  briefly  discuss  the  results  of  relevance  to  spacecraft -space- 
plasma  interactions. 

Starting  with  (1)  and  (2)  using  (3),  (4),  (5)  equation  (6)  is  obtained. 

Using  the  non-dimensional  quantities  given  in  (7)  equation  (8)  is  obtained. 
Assuming  that  for  x  <  0  the  plasma  is  not  perturbed,  while  for  x  >  0  there  is 
no  plasma  at  all  (i.e.,  vacuum)  and  assuming  that  for  x  <  0  the  plasma  is 
Maxwellian,  the  boundary  conditions  for  the  equation  (8  =  8^),  i.e.,  for; 


—  00 


+  0<3 


:r 


r(~u 2) 


^  v ,rr  ^  Jl  ■ 


p  p 

First,  the  asymptotic  cases  =—  -*•  0  and  —  ■+  »  are  discussed.  Thecaseof 
Te  Te  1  i 

—  0  (or:  6  *  —  <<  1):  This  case  is  essentially  the  'neutral  approximation' 

mirnmmmmmmmmmmk  ^rnrnmmm 

case  discussed  in  Section  II.  Hence,  expressions  (10)  and  (11)  apply  here  and 
the  distribution  function  and  n(T)  i.e.,  the  ion  concentration  distribution  are 
sought  from  (20)  with  m  0  (based  on  (25)).  Thus  for  n (C)  we  obtain: 

f  °° 

^  ^ ^  J  f  C"Ut)^U  =  jj  -  C't'jJ  (26) 

T 

Now,  fortl>>l  (note  that  large  T.  values  correspond  to  large  values  of  x) 


we  obtain: 


Note :  Expression  (27)  represents  n(t)  for  T  >>  1  (or:  T  +  “)  near  the 

vicinity  of  the  body  or  plane  behind  which  there  was  initially  a  vacuum,  and 

the  plasma  expanded  into. 

T  X 

The  Case  of  ■=—■•*  (or  6  ■  —  >>  1)  :  In  this  case  the  thermal  velocity  of  the 

mmmmmmrnmmm7 ±  wmmmm 

ions  is  less  than  the  velocity  of  the  ions  which  corresponds  to  the  ordered 
motion  (stream  velocity)  produced  by  the  electric  field.  In  this  case  (argues: 
Gurevich  et  al . ,  1969)  there  is  no  need  to  use  the  kinetic  equation.  For  this 
case  the  equation  for  the  velocity  (v)  corresponding  to  the  ordered  motion 
should  be  adequate.  Therefore  equations  (13)  apply  here  and  from  which  follow: 

wTf I 

:  u.#/  J  ” 

where:  f  =  (and  for  one  dimension  we  are  concerned  only  with  x  as  is  the 
J  Rq 

case  here) . 


Finally  we  obtain  from  (2  8)  : 


(29' 


=  0 


dH  w  oLd 

'  r 


=  0 


Again  (29)  describes  the  ion  motion  as  does  the  equation  of  motion  for  an  iso¬ 
thermal  ideal  gas  (with  temperature  T  )  in  hydrodynamics. 

The  system  of  equations  (29)  has  two  solutions.  One  which  is  trivial  i.e., 
N  =  const,  u  =  const,  and  the  second  is: 


k-c)  =$  to  -  c 
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m 


-ftrr,  X  v>1 


where  C  is  a  constant  (C  =  0.7). 


u  *  0)  and  for  X  +  00  (n  ■  0)  is: 

U=0 , 71=1  cui 

a=z:*k  =**?(-& -z-i)  4 

And  the  electric  field  appearing  as  a  result  of  the  flow  has  the  form: 


where  I 


Y-  C“C)  - 
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V5  y  -T 
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and  the  basic  definition  of  ?(x)  is: 


Va 


(30) 


The  continuous  solution,  satisfying  (2  5)  forT  +  -  »  (X  *  K^,  or  n  »  1; 


(31) 


(32) 


(33) 


It  should  be  realized  that  the  line  f  -  -  —  (see  Figure  A-3)  is  the  line  that 

/2 

separates  the  region  occupied  by  the  gas  at  rest  and  the  region  occupied  by  the 
moving  gas.  The  first  derivatives  of  u  and  N  have  a  discontinuity  along  this 
line.  A  discontinuity  in  the  electric  field  strength  means  that  there  is  a 
charged  electric  layer  along  this  line.  This  layer  is  of  the  order  of  a  Debye 
length. 

An  important  point  to  note  is  that  for  large  values  of  (t) ,  the  expression 
(31)  gives  much  higher  values  than  those  obtained  by  expression  (27).  Kamely: 


n 


rather  than: 
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-e* f  (- 
«-Yp  C 


(34) 


This  result  is  due  to  the  strong  influence  of  the  electric  field  on  the  ion 
motion  i.e.,  due  to  the  acceleration  of  part  of  the  ions  bv  the  electric  field. 
The  general  case  i.e.,  any  B; 

Mathematically,  for  the  general  case  i.e.,  (any  value  of  £)  one  seeks  the 

solution  of  equations  (8-9)  namely,  the  determination  of  the  behavior  of  the 

characteristics.  This  is  dealt  with  in  great  detail  in  Gurevich  et  al.  (1969) 

and  Gurevich  and  Pitaevskii  (1975) . 

Here  we  will  bring  only  the  results  of  interest  to  us. 

y 

Figure  A-4  shows  the  relation  between  n  *=  plotted  as  a  function  of  (t) 

/  x\/ Mv  V-P  T  0 

PC  for  various  values  of  S(r  vf)  .  The  dotted  line  is  the  case  of 

isothermal  hydrodynamics  which  corresponds  to  the  limit  S  ■*  K  (expression  (30)) 

It  is  seen  that  with  the  increase  of  £  the  distribution  [nCC)]  approaches  the 

hydrodynamic  case  given  in  expression  (30)  . 

Figure  A-5  depicts  the  same  on  a  logarithmic  scale  (for:  £  *  1).  The 

dotted  line  is  obtained  using  the  free  molecular  motion  (£  -  0)  i.e.,  the 


A 
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heutral  approximation'  (given  in  expression  (27)). 

N 

From  Figure  A-5  we  see  that  for  larger  values  of  t  the  values  of  n(*  — ) 

N 

o 

decrease  slower  than  do  the  values  of  (n)  in  the  'neutral  approximation'. 
This  is  due  to  the  influence  of  the  electric  field. 

The  asymptotic  behavior  of  n(XD  for  T  >>  1  is  given  by  (31)  namely: 


n(t:)  ~  c  f-x-fal') 


(35) 


with  C  »  0.7  for  B  “  1.  For  arbitrary  values  of  6  the  'constant'  C  is  C(8) 


and  is  given  by: 
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(36) 


Figure  A-6  shows  the  distribution  function  g  as  a  function  of  u  for  different 
values  of  Cc)  and  for  B  ■  1.  This  is  a  very  important  figure  since  it  shows 
quantitatively  how  the  electric  field  influences  the  distribution  function  of 
the  ions.  As  seen  from  this  figure,  for  X  <<:  -1  the  function  g  is  an  unperturbed 
Maxwellian.  Thus  for  example  for  T  ~  -2  the  distribution  function  g  is  close 
to  [exp(-u2) ]  .  And  in  fact  for  X  <  0  these  functions  are  close  enough  to  the 
distributions  obtained  from  the  "neutral  approximation".  This  however  is  not 
the  case  for  large  positive  values  of  (X) ,  or:  for  T  >  1.  Figure  A-6  shows 
the  gradual  conversion  of  the  Maxwellian  distribution  function  to  a  f-tvpe 
function,  i.e.,  to  a  needle-type  (or:  needle-shaped)  function  is  the  region  of 


"strong 
where:  t 


rarefaction"  i 


•  e . , 


in  the  region  where  t  >>  1  (and:  T 


Again  the  needle-shaped  functions  (for  increasing  values  of  X)  indicate 
the  strong  acceleration  of  a  portion  of  the  ions  which  fill  in  the  rarefied 


region . 

Moreover,  Gurevich  et  al .  (1966)  show  that  the  mean  energy  of  the  ions 
for  large  values  of  X  is  equal  to: 


ro 


T-u' 


(37) 


from  which  follows  that  the  energy  of  the  ions  can  exceed  the  initial  thermal 
energy  by  orders  of  magnitude.  The  effect  of  acceleration  can  even  be  larger 
for  a  bithermal  (e.g.|3  >1;  ">  T^)  plasma.  Note  that  this  latter  case  is 

common  in  laboratory  simulation  studies. 

Furthermore,  Gurevich  et  al.  (1966,  1969)  and  Gurevich  and  Pitaevskii 
(1975)  show  that  as  X  increases  the  thermal  velocity  spread  of  the  ions 
diminishes  rapidly.  And  that  the  quantity:  ^(Tj  effectivej*=  2T[u  -  u  ] 
which  has  the  meaning  of  an  "effective"  ion  temperature  and  varies  with 
increasing  T  in  the  form: 


( 


T  -ay-p  (-2,-r/a) 


(38) 


The  latter  has  an  important  effect  on  the  stability  of  the  region  behind  a  rapidly 
moving  body  e.g.  the  wake  of  an  ionospheric  satellite.  It  should  be  realized 
that  this  decrease  in  the  temperature  refers  to  the  thermal  motion  of  the  ions 
along  the  x-axis  only.  Since  the  plasma  is  collisionless,  the  temperature  of 
the  ions  in  the  flow  direction  is  not  changed  in  the  expansion  process. 
Consequently,  for  large  values  of  (t)  the  electron  temperature  may  be  several 
orders  of  magnitude  greater  than  the  effective  ion  temperature.  This  then 
reduces  the  stability  of  the  plasma  (i.e.,  since  f  T^ef  f  ective) ) . 


In  all  the  above  discussion  Z 


x  _ o 

z\  2T 


where:  =  the  velocity  of 


the  body  (which  is  supersonic) ,  z  is  the  coordinate  in  the  direction  of  motion 
of  the  body  (z  <  0),  x  is  the  coordinate  in  the  direction  of  the  normal  to  the 
edge  of  the  body  and  the  results  are  applicable  for  the  region: 
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A  detailed  analysis  of  the  stability  in  the  near  wake  zone  was  carried 
out  in  Section  12  of  Gurevich  et  al .  (1969). 

For  the  case  of  an  expansion  of  a  plasma  behind  a  plate  (or: rarefied 
plasma  flow  behind  a  plate)  the  region  of  stability  is  shown  in  Figure  A-7 . 

In  this  figure  the  results  for  a  plasma  with  6=1  are  shown.  Ion  acoustic 
waves  are  unstable  in  the  region  above  the  curve,  and  the  region  of  instability 
expands  with  the  increase  in  8.  Note  that  the  vertical  and  horizontal 
coordinates  in  this  figure  (p,  t)  are  defined  as: 


where  z  is  the  axis  parallel  to  the  velocity  of  the  incident  plasma  flow  (flux) 


and  the  x-axis  is  orthogonal  to  the  velocity  (V  )  and  to  the  axis  of  the  plate. 
The  plate  has  a  width  of  2R^  which  is  infinite  in  the  y-axis  direction. 


Flow  of  a  plasma  containing  a  mixture  of  ions 


The  case  of  a  plasma  (e.g.  the  lower  ionosphere)  having  more  than  one  ion 


species  was  dealt  with  by  Gurevich  et  al.  (1969)  and  expanded  in  Gurevich  et  al. 
(1973).  Briefly,  if  we  have  a  plasma  with  two  kinds  of  ions  (e.g.  [0+]  and  [H+] 
with  atomic  mass  units  and  y^  than  the  following  equations  should  be  solved: 


~kw,  k  [nJIM 
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The  boundary  conditions  for  "C  -*•  -  ®  are: 

<?,=  e*f  (-pO 
$*=«■  e>f(-K'R: ) 


(40) 
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The  solution  of  the  system...  depend  now  not  only  on  £  «=  —  ,  but  also  on  the 

^  N 2  o 

ratio  of  the  masses  (— -)  and  on  the  ratio  of  the  initial  concentrations  a  **[— — ) 

M2  vMo 

If  >  K.  (say;  H(0+)  =  16  and  11(R+)  ■=  1)  than  at  large  values  of  (t)  the 

lighter  particles  dominate,  irrespective  of  their  initial  concentrations.  Hence, 

if  we  assume  that  <<  (which  is  the  case  of  the  terrestrial  ionosphere 

for  H  <  (400  -  500)  km)  we  find  that,  behind  an  ionospheric  satellite  with 

M2  + 

—  =  16,  the  concentration  of  [H  ]  ions  is  much  larger  than  ttie  concentration  of 
M1 

[0  ]  even  if  the  [F.+]  ions  form  less  than  (1  -  10%)  of  the  ions  in  the  un¬ 
disturbed  plasma. 

D.  About  the  flow  over  bodies  having  a  similar  cross  section  to  the  flow 

The  question  of  the  flow  over  bodies  having  different  shapes  and  cross 
section  was  addressed  by  several  papers  of  Russian  group  (e.g.  Al'pert,  1976; 
Al’pert  et  al.,  1965;  Gurevich  et  al.,  1969).  Briefly,  they  show  that  based 
on  the  law  of  similarity  the  pictures  of  the  flow  past  bodies  are  similar  if 
the  bodies  have  the  same  contours  in  the  variables  ( p ,  t)  where:  j>  **  — 

r  *  2Tp  w  z  i  v-  ° 

and  t  »  (—-)  *  ( - — )  7-  (p-)^  where  • 
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thermal  ion  velocity. 
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Note:  The  laws  of  similarity  found  for  the  rarefied  plasma  flows  (e.g.  Gurevich 
et  al.,  1969)  and  analogous  to  the  laws  of  similarity  is  supersonic  aerodynamics . 


In  fact  in  dealing  with  the  case  of  rarefied  plasma  flow  over  a  disc  of 


radius  Rq  it  is  argued  that  the  disc  case  is  of  general  significance  since  it 
represents  the  flow  over  any  body  of  rotation  provided  that  the  length  of  the 
body  (L)  is  not  too  large,  i.e.. 
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Is. 


<?<  R, 


(42) 


and  that  the  results  obtained  for  the  disc  will  be  valid  for  such  a  body  at 
sufficiently  large  distances.  This  is  particularly  so  for  flow  past  a  body  of 
rotation  with  a  flat  near  side. 

O  Comments  on  the  Influence  of  the  electric  field 

The  influence  of  the  electric  field  (Gurevich  et  al.,  1969;  Al'pert,  1976) 
is  discussed  via  a  comparison  with  the  'neutral  approximation'  (where  no  elec¬ 
tric  field  effects  exist). 

Note :  In  Al'pert  et  al .  (1965)  the  expressions  used  for  the  ion  densities 
are  in  fact  those  of  the  "neutral  approximations". 

For  example  for  the  case  of  flow  over  a  semi-plane  the  difference  between  the 
"neutral  approximation"  and  the  calculations  of  ion  distribution  (based  on  the 
discussion  given  in  this  appendix)  is  shown  in  Figure  A-5.  As  seen  from  the 
figure,  the  distributions  of  ions  and  neutrals  are  similar  for  negative  and 
small  positive  values  of  T  * yP/t  <  0.5,  in  fact  they  are  identical  to  within 
10%.  On  the  other  hand  in  the  region  of  maximum  rarefaction  (p  =  0  or  large  T) 
they  differ  by  orders  of  magnitude.  Gurevich  et  al.  (1969)  have  performed 
similar  comparisons  for  wakes  behind  a  plate  and  a  cylinder  and  reached  similar 
conclusions.  Furthermore,  for  the  case  of  a  sphere  (or  a  disc)  the  difference 
between  ions  and  neutral  particles  in  regions  which  are  not  strongly  rarefied 
(e.g.  N  >_  0.2IO  the  deviation  of  the  distribution  of  ions  from  the  distribution 
of  neutrals  are  more  significant  (see  Figure  A-8) .  The  role  of  the  electric 
field  increases  substantially  with  the  increase  of  the  ratio  [T^/T^]  (see  Fig. 

i 
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A-8,  for  8  «  4) .  Another  example  showing  the  difference  the  distribution 
according  to  the  "neutral  approximation"  (or:  free  motion  of  molecules) 
and  ions  in  the  wake  of  a  plate  and  a  cylinder  is  given  in  Figure  A-9 . 

•  About  theory-experiment  comparisons  and  the  use  of  the  expression 
[N/N^]  (given  in  the  text) 

Gurevich  et  al.  (1969);  Gurevich  et  al.  (1973)  and  Al'pert  (1976)  have  6hcvn 
results  for  theory-experiment  comparisons.  The  first  two  references  performed 
comparisons  with  in-situ  measurements  (mainly  using  results  by  Samir  et  al . ) 
and  Al'pert  (1976)  showed  comparisons  with  both  in-situ  and  laboratory  measure¬ 
ments  . 

In  this  section  we  will  show  the  degree  of  agreement  obtained  via  presenting 
a  set  of  figures . 

* 

Figure  A-10  shows  the  degree  of  agreement  between  measurements  of  I  •  f(g) 
obtained  from  the  Ariel  I  (U.K.)  satellite  and  computation  by  Gurevich  et  al . 
(1969)  using  an  expression  for  flow  behind  a  cylinder. 

This  expression  is  formula  (54)  p.  830  in  the  above  reference.  The  compu¬ 
tation  was  done  for  8=1  and  b  *  3.75,  and  is  displayed  by  the  solid  line.  As 
seen  there  is  good  agreement  for  8  <  120°  only. 

Figure  A-ll  shows  the  degree  of  agreement  between  I  £  =  f(8)  at  a  distance 

of  5R  from  the  center  of  the  Ariel  I  satellite,  as  measured  bv  the  boom  mounted 
o 

electron  probe  (guarded  plane  probe)  and  between  the  Gurevich  et  al .  (1969), 

Al'pert  (1976)  computations  given  by  the  solid  line.  Here  as  in  Figure  A-10 

the  computation  was  done  for  £  =  1;  b  =  3.75  using  the  expression  for  the 

"neutral  approximation".  Although  the  measurements  were  for  the  current  ratios 

I*e  =  [I  /I  J  (which  is  the  ordinate  of  the  figures)  anc  the  computations  for 

N*  =  [-r~]  it  was  assumed  that:  [I  /I  ]  s  [N  /V.  ].  Of  particular  interest  in 
b  e  o  e  o 

o 

Figure  A-ll  is  the  fact  that  the  "neutral  approximations"  (at:  Z  =  5?.^  from 
the  center  of  the  satellite)  predicts  satisfactorily  the  degree  of  electron 
current  depletion  at  the  maximum  rarefaction  ( ^  =  0;  large  TO  on  the  wake  axis. 


J*  t- 
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Figure  A-12  shows  the  comparison  between  I  ■  f (0 )  measured  by  an  elec¬ 
tron  probe  flush  mounted  on  the  surface  of  the  Explorer  31  satellite  and  the 
Gurevich  et  al.  (1969)  computations. 

Note:  The  probe  on  the  Explorer  31  was  similar  to  that  on  the  Ariel  I  satellite. 

The  result  shown  is  similar  to  that  shown  in  Figure  A-10.  An  the  expression 
used  in  the  computation  is  identical  to  that  used  for  Figure  A-10. 
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Figure  A-13  shows  a  comparison  of  some  Explorer  31  results  with  results  of 
computations  using  expression  (1)  given  in  text  of  the  paper  i.e.,  the  expression: 
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As  seen  from  this  figure  the  agreement  between  theory  and  experiment  is  very 
reasonable. 

If  [H+]  is  the  major  ionic  constituent  for  altitudes  above  (500  -  600  km) 
in  the  terrestrial  ionosphere  than  the  effect  of  the  electric  field  on  the  motion 
of  the  ions  is  small  (since  the  ionic  Mach  number  S(H+)  is  small)  and  in  the 
boundary  of  the  wake  the  field  can  be  ignored.  The  angle  (y^)  in  expression 
(1)  characterizes  the  wake  boundary  as  seen  by  the  probe  situated  at  on  the 
satellite.  For  the  comparisons  shown  in  Figure  A-13,  the  angle ■=  45°. 

Another  conclusion  (or:  recommendation)  of  Gurevich  et  al.  (1969)  is  that 
when  the  relative  value  of  [H+]  in  the  plasma  exceeds  about  30%,  expression  (1) 
can  be  applied. 

j  /  jQ.  _  igo°J  4. 

Figure  A-14  shows  the  dependence  of  [- — ^j0 - "]  on  n^+I  N  (H  )  /K^total)  ] 

in  the  range  0  <C  n^+  <,0.75  (after  Gurevich  et  al.,  1969),  in  the  quasi-neutral 

region.  The  solid  line  shows  the  results  of  the  calculation  using  expression  (1) 

for  T  «  constant,  the  dotted  line  shows  is  essentially  the  same  but  taking 
e 

account  of  the  variation  in  Tg  and  the  points  are  Explorer  31  data  points  (after 
Samir  and  Wrenn,  1969) . 

In  Gurevich  et  al.  (1973)  it  was  shown  that  in  a  plasma  which  expands  into 

a  vacuum  a  portion  of  the  ions  is  accelerated  and  the  energy  they  may  obtain  by 

X  V 

the  process  of  acceleration  can  be  of  the  order  of  (10  -  10  )  T  .  In  that 

publication  the  excitation  of  waves  (e.g.  ion-acoustic  waves)  is  also  investigated. 
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The  treatment  is  that  of  a  plasma  having  two  ionic  species  i.e.,  [H+]  and  [0+]. 
Based  on  the  computation  given  in  this  paper  the  behavior  of  Figure  A-14  was 
replotted  and  as  can  be  seen  in  Figure  A-15  there  is  a  better  agreement  between 
theory  and  experiment.  In  Figure  A-15  the  dotted  line  represents  the  calculation 
using  the  expression: 


i.e.,  essentially  using  expression  (1).  As  is  seen  from  this  figure  the  solid 
line  depicts  a  better  agreement  with  experiment  for  n^.  >  0.4. 

In  conclusion:  It  appears  that  for  specific  cases  of  practical  interest 
in  the  terrestrial  ionosphere  and  in  planetary  ionospheres  semi-analytical 
expressions  for  the  distribution  of  particles  around  a  body  may  be  adequate 
and  there  may  be  no  practical  need  to  seek  very  elaborate  and  costly  numerical 


codes . 

This  conclusion  may  also  be  of  practical  interest  to  spacecraft-charging. 


Fig.  A-!  .  Classification  of  initial  discontinuities  Different 
relations  between  the  bound.tr>  distribution  iuiudont  y 
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Fig.  A-2.  Non-dimensional  force  F(r)  for  the  same  cases  as 
in  big  I.  The  sign  of  the  velocity  t  is  different  in  the 
curves  I  and  2  of  casetb)  Curve  2  of  case  (c)  corresponds 
to  the  fuc  expansion  of  plasma  into  vacuum:  N,  -  0 


Fig  .  A-4  .  Flow  past  a  semi-plane  The  dependence  of  concentration 
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for  different  values  of  /)  =  7V/7V.  The  dotted  line  is  obtained  in  the  case  of  isothermal  hydrodynamics, 
which  corresponds  to  the  limit  /?-*  co. 


Fig  .  A-5  .  xhe  same  on  a  logarithmic  scale  (fi  —  I).  The  dotted  line  is  obtained  for  free  motion  of 
molecules,  which  corresponds  to  the  limit  0-+O  (’neutral  approximation'.) 


Fig  •  A-6  .  Flow  past  a  semi-plane.  The  distribution  function  g(u)  for  different  values  of  r. 


F ig  .  A-7  .  The  boundary  of  the  region  of  instability  for  flow  past  a  plate  (0  -  I)  Ion  acoustic  waves 


are  unstable  in  the  region  above  the  curve. 


